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THE principal object of the calculations contained in the following paper, is to 
investigate the figure which a fluid, consisting of portions varying in density accord- 
ing to any given law, would assume, when every particle is acted upon by the attrac- 
tion of every other and by a centrifugal force arising from rotatory motion. To 
what extent this may have been the original condition of the earth, is a doubtful 
question ; and although observation does not fully warrant this supposition of the 
regular arrangement of parts, it has necessarily been made the foundation of most of 
the mathematical calculations connected with the investigation. Before proceeding 
to this problem, it is necessary to calculate the attraction of a body of any given 
figure, and consisting of strata, varying in their densities according to any given law; 
and it is in this problem that the principal difficulty lies. The elegant method of so- 
lution discovered by Laplace is well known ; and I have followed his steps as far as 
the point where the equation, known by his name, first appears. In order to illus- 
trate the nature of the deviation which I have there made, it will be necessary to 
mention some of the principal steps of the two methods. 

By means of a theorem, which Laplace laid down as true of all spheroids that 
differ but little from spheres, and the properties of the integral of the equation re- 
ferred to, he was enabled to substitute the easy rules of differentiation for the more 
complicated inverse processes, and thus to compute the attraction of that class of 
figures. It has, however, been since discovered by Mr. Ivory, that this theorem is 
true only of spheroids of a particular kind ; and, consequently, to this kind the solu- 
tion of the problem is restricted. This defect, and the indirectness of his analysis, 
led other mathematicians to consider the question; and, in 1811, Mr. Ivory pub- 
lished his method, which has the great advantage of being more direct, though 
equally limited. 

The method given in the following paper does not appear to be confined in its 
operation to any particular class of spheroids ; since the coefficients of the series, into 
which the required function is developed, are determined absolutely, without any re- 
ference to the form of the spheroid to which they are about to be applied. The prin- 
cipal change consists in the different manner of treating this partial differential equa- 
tion. Laplace and the subsequent writers on this equation, both as applied to the 
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calculation of attractions and the mathematical theory of electricity, suppose the co- 
efficients of every term of the series to be expanded into another series of the sines 
and cosines of multiple arcs ; and they avail themselves of the property which these 
terms possess of vanishing, in certain cases, when integrated between certain limits. 
The success of this plan, however, depends upon the restricting hypothesis above re- 
ferred to, that the radius vector of the surface of the body is capable of expansion in 
a series of terms, each of which satisfies Laplace's equation. The following method 
shows that the coefficient of the general term of the first series is independent of one 
of the variables, and thus dispenses with the second series of expansions. This re- 
sult I have arrived at, by first obtaining the integral of Laplace's equation in its 
most generalform, and deducing the arbitrary functions introduced therein, from 
considerations which enter previous to the equation of the surface of the attracting 
body. These coefficients being known, it is evident that the attraction of any homo- 
geneous body on a point within or without it may be immediately found when the 
equation of its surface is given, since it then depends only on a series of explicit and 
definite integrations of known functions, which can always be effected, at least ap- 
proximately. From this, the attraction of a heterogeneous body, similarly circum- 
stanced, may be found by the usual method of dividing it into concentric layers, and 
summing the several attractions of these, deduced as above. 

By Substituting the attraction so obtained, in the equation of equilibrium of a fluid 
body, Clairaut's theorem is immediately deduced; and, from a peculiarity in the 
functions representing the attraction, it will be seen, that the same principles with 
longer processes may be carried on indefinitely, without the necessity of actually 
determining the precise form of those functions. 

The restricted species of spheroid above referred to, comprises all surfaces of revo- 
lution ; so that it is sufficiently extensive for most practical purposes ; but the inte- 
gration of Laplace's equation renders the analysis more direct, and the theory more 
complete. 

On the General Problem of Attractions. 

1. Let g represent the density of a body at the point (#, y, z) ; and let /, g, h be 
the coordinates of a particle attracted by the body, parallel respectively to the axes 
x,y,z\ then, if the power of attraction be inversely as the square of the distance, the 
resolved part of the attraction of the body, parallel to 

x is ffC q{f-*)dxdydz 

JJJ {(/-*) i + (y-y)* + (A- *)*}*' 



. Crr g(ff — y)dxdydz 



{(/-*)■ + fo-y)* + (A-*) 8 }- 
™fff 



g(h — ss) dxdy dz 



{(/-*) , + («'-y) 9 + (*-*) 8 }*' 
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the liiftits of integration being determined by the equation to the surface of the 
body. 

2. Let V represent the sum of the products of each particle by the reciprocal of its 
distance from the attracted point. 

Then V -JJJ { (/, f )» + ^yy + .( i ,.^ ; and > b y differentiating V, we obtain the 

well-known property -jj* + ~&tf + ~Jjfi =: °> or — 4 >r §', according as the attracted 

particle is not or is within the attracting mass ; g f being in the latter case the density 
of the attracted particle*. 1 By transforming these equations to polar coordinates, we 

obtain 

d 2 V ■ 2 dV , 1 d*V , 9 _dV ■ 1 d 3 V . ■ , 

^^ i — *— -l. — _-— j. r^cot# -J-* + - o - 2a Tir = 0, or — 4tp, 

dr 2 ' r dr x r 2 dv 2 ' at/ * r 2 siir0 d<p y » ' 

and 

""l/o */ o */o {r 8 +r f8 — 2rr' (cos cos 0'+ sin sin 0' cos (<p — $ / ))}*' 

^ a 

where r 2 — / 2 + g 2 + & 2 , cos = ,, ~ 2 pr , tan <p = y ; and similar expressions 

are true of r' 6 1 and <p- in terms of x p y> z. 

Put cos 6 = p, and cos 6 l = ^ y and they become 

r^ + ^U 1 -^ 2 )^)*!^!^ ( L > 

/*** /*1 /»£*• g^dr* dp!dtf 

V —Jo J -I Jo { r *. + r ^ _ 2/V (jxjx' + a/(1 - j* 2 ) V(l -V 2 ) cos (<p — <p'))K * ^ 

3. Expansion by the binomial theorem shows that 

{ r 2 4. r '2 _ 2 r r' (^ r + V(l — f* 2 ) 1/ (1 — fi/ 2 ) cos (<p - p'))}* 

may be expressed either in powers of r or ofV-; thus 



.rc i ^f y/2 ^,m 



where P n is a symmetrical function of p, V(l — /a 2 ) cos <p, v"(l — ^ 2 ) sin <p on the one 
part, and y! 9 i/(l — V 2 ) cos p', •(! — ^' 2 ) sin 0' on the other. 

By substituting the first expansion in (2.), and the value of V so obtained in (1.), 
we have a series of equations 

d * V n , / , nB 1 d? J d[*.'d<p ! 



fXjr t {£(<' -*>£ ) + r^ + «(» + dp.} 



^7 «~l 



= 0, or — 4^^, 
except when n = 2 ; and in all cases 

which is the equation of Laplace's coefficients. 

* Vide Pratt. Mec. Phil., § 168. Laplace, Mec. Cel. liv. iii. t Vide Pratt. Mec. Phil., § 169. 
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4* This equation was not integrated; but by a skilful use of its properties, 'the' 
problem of attractions was greatly simplified by Laplace. He laid down a theorem, 
respecting the surfaces of all spheroids of small deviation, that their radii vectores 
might be developed into series, every term of which would satisfy the above equa- 
tion ; and he also gave a method of expansion. By means of this theorem, the pro- 
blem could be solved for spheroidal bodies which differ but little from spheres ; but 
its generality has been greatly restricted by the researches of subsequent writers *, by 
whom it has been shown that it is true only for bodies whose radii are expressible in 
rational and integral functions of (t) r , V"(l — • f*/ 2 ) cos <p\ \/(l — pi 2 ) sin iff. Among these 
are the ellipsoid and elliptical spheroid, and a large class of other spheroids. In these 
papers I have adopted a different proceeding ; I integrate the equation itself generally, 
and determine the arbitrary functions contained in the integrals by the circumstances 
of the problem itself. In consequence of the peculiar form which T? n then takes, V 
may be found by effecting the operations indicated, which are only explicit integra- 
tions. 

5. I shall now proceed to integrate this equation. 

Consider ^ and <p as functions of two new variables X and Y, to be determined 
from the equations, 

dY = ^(d<p + k'dt*), 
where k and k' are the roots of the equation (1 — p 2 ) k 2 + , _ a = 0. These roots 






are -=; — ~ r - } whence we obtain 



1 +£ „„a v _ *. i / — ri^^f 



X = ^ + Jx/-llog^,andY = 9-iv'-llogf^. . . (4.) 



d^ — \dpj dX 2 "•" * dp dp dXdY "*" \dp) dY* "+" dp* dX + dp 3 dY } 

_ 1 ^P. 2 d»P. _ 1 g P 2 p V-=T (d P. dP\ 

- ~ (1 - p*f dX* "+■ (l _ pZf dXdY (1 - ft 2 ) 2 dY 2 + (1 - ft 2 ) 2 \dX ~~ ~dY) 



dP n dXdP n 



,dYdP n _ V -I (dP n _dP\ 
dp — dp dX + dp dY ~~ l-p* \dX dY J ; 

rf 2 P M (dxy<pv dxdY d*p n aeywp n ,d*xdv d»Y<n>. 

df ~ \df/ tfX 2 "T z d<p d<p dXdY~T~ \d<?) dY* "^ df dX T df dY 



~ dX* "T"^ dXdY~r rfY 2 
Substituting these in (3.), we obtain 



4 d?P 

'* See two articles by Mr. Ivory in the Philosophical Transactions, 1812. 
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From (4.), by subtraction. 



— 9 



whence 

_ s -(x-Y)^^i_ 1 

and 

_ 4 s- ( x " y ) ^~ >• 



I 2 ' — ______ ... 

i - fb - ^ (x _ y) ^—j + ^ _ ^ I (x z Y) v^T + gi(x I Y ) v^i) 2 — oX- Y' 



cos 2 



consequently 

d 2 P w . n(n -f 1)P„ ^ T ' t x 

^^Y+7 2 X--Y = » Letn(n+1)F^ 

4 cos 2 - 



Let 



or 



Let 



Then 



and 



2 



^ = ^ then^ + x- : Y = °> andP » = - «7 cos2 T"' 

2 

<2P % d*v 4 9 X~Y <2-t? 4 X~Y . X-Y 

v =dY= ~-dTdY^ C0S ~2 dX ¥ cos -2~ sm ~2~ ' 

d 2 v dv X—Y av 

dXdY + dX tan ~2~ + " TX^Y = °- 

4 cos - ' 



dt>. . 4 X-Y 

5Y + »tan.- ¥ - = / i . 

d 2 t> , dv X-Y g d* 

dXdY + dX tan ~~ T~ + " ^=Y - __.' 

_2 cos ■ 

2 

dt a-2 o ■■ 4 _X-Y'_* . 

rf_ + -4 ;X-Y = ' or i> = - ^Ts cos2 "^-Jx ' 



cos 2 



2 



whence 



rfw 4 «X— Y (Ft 4 X-Y . X-Y_* 

dY = - ^^2 COS 2 -^— ^rfy ~ JIT. COS -g— Sin -y- -^ , 

4 -X-Y d?t 8 X-Y . X-Yrf* 

' = -^T2 C0S 2 — F" JX^Y ~ «^2 cos -T- Sm — — dX' 

or 

_ft_,^ . X-Y , a-2 

dXdY "+" rfX J tan "2 f" * " 2 X-Y — ".°- 

4 cos^ 

2 

Let 

d* , „ _,. X-Y 
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and by repeating a similar process, we obtain 

dXdY~r dX 6l3,n ~ h # 2 X-Y~ U * 

4 cos 2 

2 

By observing the assumptions here made, and the results obtained, we find that in 

(dV \ X — Y 

the first assumption \^y ~ v )> ^ e coe ffi c i ent of P^ tan — ^ — is 0; in the second, 

2£ y 

that of v tan — ^ — — * ; an d so on, in the order of the natural numbers ; and in the 

ft ft _— Q 12 •""-" fi flJ ■"■""• /2 (ft mmm 1 r 

results, the numerical coefficients are 1, -^; 2, 4 ; 3, 4 . .. generally rc, ~ -. 

I shall prove this in the general case, by showing that if it is true of one value of n 
(as we see it is), it is true of the next value, and so on. Let the (n — l)th substitu- 
tion give 

(Pg .dg f . X-Y , a-(w-l)(w--2) _ A 

XdY + dX^" ! ^ tan 2 + £ . J-Y "^ U5 



4 COS 3 



and let 



-x+ ( w _ l) ? tan— — = *: 



then, as before, 

d 2 g , v X— Y # g w — 1 1 __ rf s . 

<ZX^Y + ^~~ IJtan— g— '^x + — 2~ 2 X-Y^-rfX 9 

cos 2 

At 

and, therefore, 

ds . a — (n — \)n ^ t ds 4 X — Y 

* t3v = 0, and g = - ^ ^/-iw CQg2 

4cos a 



rfx + r 2 x-y = °> and e = -rfX a-( w -i) W C0S ' — 



XV Y V 

— ~ X . .A. — X 

'■«- ,r , 4C0S « S111 ^ 

<2<? d?s 4 „X-Y <?s 2 2 



COS 2 



dY _ dXdY <*-(«- 1)^" 3 2 tfX a-(n-l)n 

Consequently 

<Fs 4 ,X-Y ds 4 X-Y . X-Y 

cos 2 • — 7T- tv 7 rtr cos — z — sin — ^ — 



or 



dXdYa-(n-l)n^° 2 dX a-{n-l)?i wa 2 

ds 4 (a- 1) X-Y . X-Y _ 

~ dXo-{ti-l)n C08 2 S,n 2 "■*» 

d*» 4 „X-Y , ds 4» X-Y . X-Y ■ 

cos 2 ■ - . + -pf - — t-— T V- cos a sm — 5— + * = 



rfXrfYa-(»-l)ii 4 '" a 2 -T dXa-(n-l)ti 2 

that is, 

d 2 s , ds x X-Y ■ a-n(n-\) „X-Y „ 
iX7Y + ^X wtan ~2-+*-4 COS 2 -^-^^ 

and, therefore, the law of coefficients, as above stated, is correct. 
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Restoring the value of a, we get 



d?s ds X — Y , s.n 1 

dXdY + dX ntan 2 +1 7T^Y~°> 



cos 2 



ds X — Y 

which is integrable ; and its integral is ^y + n • s tan — § — = some arbitrary func- 
tion of Y, as % Y. Integrate again, then 



S = Z J 2 \J * J 2 ^YdY + ^Xj, 

where ^ is arbitrary. Effecting these integrations, and reducing, 

— 2»X— -x / / 2wA-— I xr 7 _ r \ 

S = C0S — 2~\y c0S ~ 2~-%Y<£Y-f^X). 

To return to P B , we have the following systems of equations : 
P n =fvdY, 



t? = g 



£ = COS * — 2~~y? C0S 2~~'^ Y 



-2(w-i)Y — X /» 2 ; (» - 1) Y — X JAr 
g= cos — g— y ,9 COS v ' — g— .tfY, 

*= cos- 2w ^^(ycos 2w ^^. % YrfY + 4x); 

whence 

Now cos 1=5 = cos (/Tl log \Zt^J) = J (\Zt="? + a/tt£)> 

Y — X 1 

and cos 2 — ^— = L — - 9 ; and the complete integral will be expressed, by substituting 

for X and Y in terms of ^ and <p. 

6. But an important point yet remains to be determined. The original equation, 
being a partial differential equation of the second order, can only involve in its in- 
tegral two arbitrary functions. But here, after % Y and 4 X have come in by two inte- 
grations, we have n integrations to perform with respect to Y. It would seem, there- 
fore, that no constant or arbitrary function of X must be added in these integrations. 

MDCCCXLI. M 
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Such is not the case. At each integration a function of X must be added, and thsse 
functions determined by reference to the original differential equation*. 
7. Returning to the value of P % , we have 

p. = * (» -{f=-i log m ) + + (» + 1 ••=! % ;-±#. 

Now in the calculation of attractions, where P is the coefficient of r° in the ex- 
pansion of 

{r 2 + r' 2 -rr f (pp! + 4/ (1 — ^ 2 ) ^ (1 - p/ 2 ) cos (?) — <p')) } ~" ^ 
we know that it is 1 ; consequently 

* (»-£•=-! log Jie). + + (* + 1 v/^T log £*) = 1, 



and expanding by Taylor's theorem, we get 

2 \2 10 & 1 - «/ 



XlP-x'i?>ix/-llog^ / * 



A* 



+ "SI */~' (i"* "r^) 3 + «*■ 



+ ^ <P + ^ P 2" V — 1 log j 



1 + /* 



/* 



- *? (y »* m 



2.3 



V—l i} log [±*)" + 



&C. 



j 



1. 



By equating the coefficients of the same powers of-g log j-37^, we have % x <p + ^ <p 

= 1 , and ^' ^ — %'x ^ = 0, or ^ p — xi P = constant. 
Therefore ^ p and ^ <p are absolute constants, and their sum is 1 ; whence it fol- 



lows that x P = 0. Let %// <£ 



2> 



then 



w 



V • • • v 



2Y-X 



COS — 2~~J C0S " r ~ * ~ C0S 



2 Y-X/n 

2 




9 y — X 

— 2 — - axax...(n times 



>) 



Effecting these integrations, we find that P n consists of a series of powers of 



* The common differential equation (1 -— ju 2 ) 



n 



LetP 



n 



2 jj — -^ + n {n + 1) P« = will illustrate this point. 
-, and after substitution, differentiate n times ; then (1 — jw 2 ) — - — 2 (ft -f- 1) p - T — = 0, whence 



cfiP, 

<^ 3 



djjb-n 



z 



Jc 



dp* 



d 



P 



n* // #/ 

1___ -f- w. It is clear that no more arbitrary constants than & and m can be introduced ; and 

(l-.ps)»+l- y 

yet if the integrals were left indefinite, we might obtain an integral of an expression which should differ from 
the integral of the same expression obtained by a slightly different process, by a constant. By another inte- 
gration this would cease to be a constant, and we should obtain thus different values for P w . The fact is, that 
constants must be added at each integration, and recourse had to the original equation, to determine them in 
terms of m> k, and .p. 
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o 



tan — g— , whose coefficients may, for anything we yet know, be functions of X. 

The following process shows that constants as coefficients will satisfy the original 
equation, and determines them. The integration itself gives the coefficient of the 
highest power. 
Let then* 



tan" ±=A + d tan» ~ 1 ^=^ + c" tan* ~ 2 ^r-^ 



P =c 

n 



W 



2 



+ . • • + c (ra - 2) tan 2 ^^ + c^- 1 * tan ^-^ + c« 



J 



Then 



dV 



n 



i 



dY 



m — 2 



[ft-1] 



tan 



w — i 



Y-X . »-l 



wImm 



(/tan 



n -2 



Y-X 



L=X +....+ I c (— 2).^ y-x 



= c i 



+ —-? c" tan" " 3 



2»~ 2 .xi Y-X , »-l ,,^ J-X 
tan" + * — s 1 5 — c tan" 



(l + 



tan 5 



Y-X 



) 



[II -1] 



4" \ 



2 
ft — 3 



"T I 



J// 



c"' + 



ft 



lc ') 



tan 



2 » \ 2 

Y — X 



ft — 2 „ , 2 W ~' 2 



[ft-1] 



J tan w 



Y-X 



tt-2 



2 



MfH» 1 



»ri c » + !Lz_V)tan-» 



") 



Y-X 



+ ... + (§ c (— *>+!■ c(»- 5 ))tan^ 



Y-X 



+ (t c(W " 2) + | c( " ~ 4) ) tan3 



— X 



+ (1 c (» - 1) + 1 c (» - 3)) tan 2 



Y-X 2 

2 +2 



+4 c(» " 2 ) tan ^— ^ +.-5- c<" - x ) 



r (ra+ 1)2 



»-3 



<PP. 



W 



dX^Y 



cos" 



Y — .a. i 



n 



^— tan 

[n - 1] 2 



Y — X ft (ft -|- 1) . tt _ j Y — X 






c' tan' 



2 



-c 



m — 2 



+•=1 (•=««. +^1)^-. 



— |" m 



2 \ 2 
ft — 2 /ft — 3 



c 



C "~f" 



[n-l] 



tan 



% -3 



+ . . . + 4 (4 c(w ~ 3) +4 c(w ~ 5) ) tan3 



Y-X 

Y-X 

2 

Y-X 



+f(4 c(n - 2) +4 c(w - 4) ) 



tan 5 



XT XT 



+ 4(l c( """ 1) + 4 c( *~ 3) ) tan^=^-~ 11 c(-- 2 ) 



2 2 



* [*] = 1 . 2 . 3 . 4 . . . . w. 
M 2 
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n + l n 

n — - A — F 



4 'n 3 Y-X 

cos 3 



COS 



9 



Y x 



(n + l)2*-» ,,IJ , »(» + !) ., 



0-1] 



tan 



*wAv 



er tan 



I » * » » 



»-i 



Y-X 



H — - 4 ; c" tan w ~ ^ — g — 



4 



V. 



Whence we obtain c' = 0, c"' = 0, c v = 0, &e. 
Also 



( 



»(»+!) (»—!)(«— 2) 



)c" = 



>w — 3 



[»-2] 



.(' 



n(»+l) (k— 3)(n— 4)\ iv _ (ft — 3) (ft — 2) „ 



! 



c ; 



(?Mli) _ (»-5)(»-6) \ cVi = (n-5)(»-4) civ &c> &c> 
V 4 4 / 4 



Consequently 



P = c 



m — i 



+ 



tan" ^r^ + 



>» — "1 



2 ■ 2[ft-2](2ft-l) 
2 n- 1 



tan 



w-2 



Y-X 



2.4. [ft - 4] (2 * - 1) (2ft - 3) 



tan 



% — 4 



2 



!•• 



I"- 1 tan tt - 6 



Y-X 



+ 2.4.6.[ft-6].(2ft-l)(2ft-3).(2ft-5)~ + -- 



Now 



Y -X 



2 

2 



s/- l l °Sr~= tan- 1 (-f*</- i),... 



tan 



Y-X 



(-f«V- i); 



whence, finally, 



\ fftl 9.Tm — Sl./^W.— "H ' 9, A Vn—Al /0«_l\ /0*,_a\ "•"•••Jj 



X 



[»] ' 2.|>-2].(2ft-l) ' 2.4.[ft-4](2ft-l)(2ft-3) 

a remarkable result, showing that in this instance P n is independent of <p. 

P n being free from <p, is a perfectly symmetrical function of ^ and ^ ; and p! is a 
constant with respect to p ; therefore 

r w — ^nV [ w ] T 2 [ w _ 2 ] (2ft-l) T 2,4.[n-4] (2» — 1) (2ft-3)" r#, V 

V [«] ^ 2.[ft-2](2ft-l)^' , V* 

To determine K n for any particular value of n, we refer to the expression from 
which the two series were deduced ; namely, 

| r 2 +r i2 _ 2 rr' (^ + V(l — ^ 2 ) */(l — ^ /2 ) cos (p — ^)) } ~ ^* 
When ^ and frf are each 1, then P n = 1, which gives an equation to find K n . 
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8. Returning to the theory of attractions, we have, when the particle is internal, 

for that portion of the body which is comprised between a sphere of radius r, and the 
surface of the body, supposed to be a surface of revolution round the axis of % ; for, 
in that case R, the value of r f at the surface, is independent of <p'. 

9. Suppose, for example, we wish to find the attraction of a homogeneous spheroid, 
on a point within it. In this case g is constant, and 

a being the semi-major, and c the semi-minor axis. 

First, all the even terms vanish ; for the general even term is 

Now P 2n+1 consists of odd powers of ft/ ; and y a% + ^J 2 can be expanded 
in even powers of p! ; therefore the integral of the product (which is an odd function), 

taken from pJ = — 1 to pi = 1, is 0. AteoJ _ x P 2n +i ^P* ~ ° *• ^ t ^ ie 0( *d terms 
above the third vanish ; for the (2 n + 1 )th term is 

»-^y>..^"^=-^/-.p,.'-(4£+0""V 

-| — 2i£L /_]_ V 2n r 2 d p! 9 for it may be shown that /_ x / P t P^ dp 1 d(p f = 0, if # and/ 

be different integers. Now when rc is greater than 1, ^— ^ h -^j is a rational 

and entire function of pi, and, therefore, capable of being expressed in a series of 
Laplace's coefficients -f , the highest of which will be of the (2w- 2)th order ; and 
therefore no term of this expansion can be of the same order as P 2 n ; and the inte- 
gral of the product of any two of different orders, between these limits, vanishes. So 
the second member of this vanishes. 

* See Pratt. Mec. Phil., § 180. 

t See Pratt. Mec. Phil., § 176. Poisson, Theorie Math, de la ChaL, chap. viii. Laplace, Mec. Cel. liv. iii. 

chap. ii. 
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The first term is 



- 1 



a 2 H^- — r 2 L (e being the eccentricity). 

The third term is 

f K 2 */ Y_\ (1 - <* 2 ) (I - ^' 2 ) F ^ ^ 



* e V (y - p 2 ) ^ ( 



c« 



2 



3 tan 

(a 2 - c 2 )^ * 



3 /„« 



- l ^! "~ c ^ _ J? ™ 



3 a 2 — c 2 9 



) 



-s*e(T 



w 



! ) - 2 (t ■ 









,2 *T~ ^3 



sin- * e), for K 2 = 9. 



4 9T D f 2 ^ 

The value of V for the sphere of radius r, calculated by the usual method, is — ~ — 



consequently, for the whole ellipsoid, the value of V is 

-f--~ +2<?rg \cfi — -j— ^sm- l e - r 2 j — 3^y- ^ 2 jr 2 V3---7H ^ — sin- 1 */. 



& 



Differentiate to/, by means of the equations r 2 =/ 2 + g 2 + A 2 and ^ = — , and we 
have 



so 



so 



dV 

if 

dV 

d g 
dv 

dh 



= attraction in a? = + 2 *rg/( 1 — t H ^ — sin- x e) ; 

= attraction in y = + 2 t g> g 11 5- + e 3 ■ sin- l e) ; 

= attraction in y = + 4 sr g> A (4* ^a — e 3 — sin- l e ) ; 



which are the 
common ex- 
pressions 
otherwise 
found -j~. 



Also 
— — 

dr 



= attraction to centre =+ 4 a* g>< y+ vHif ~~ ^ /( "3"""*"^ + ~~^ — -sin- 1 e) L 

10. By a similar process, I have deduced the attraction of an oblate spheroid, on 
a point within it ; the density varying inversely as the distance from the centre. The 
corresponding expressions are 



dV 
df 

dV 
dg 



J 

+ J i ft p — 



3^_ 

4 a{ 



+ 2ff^- 



9 3 *%9 



1 - e 2 ) J e 2 ~ 3 \ 

3 V 



-e 2 ) 2 2_ 1 -e 
2e s """3 e " 



4 a (1 - e 2 ) ) e 5 



(l-<? 2 ) 2 2 1-e 
2e 3 "■ 3 e 



)log^}; 

-; Jog 1—4, 



* See Pratt. Mec, Phil., § 172. 



t Ibid. § 158, 
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^1 - m o A . 3 »8* | 1 1 / (l-gy . 2 1 - g* \ . \+e\ 

™ 77 - + 2 *<? + t^"(i-^) {? - t - vit" + y-T~7 lo snri J . 

11. When the particle attracted is external, the series in (3.) does not give a finite 
expression. Instead of taking it separately, I make it a case of a general theorem 
which follows. 

12. To find the attraction of a spheroid on a point within it, the density being any 
function of the distance from the centre, arid the eccentricity being small 

Let g <p r ! represent the law of density ; then the value of V, for the portion com- 
prised between the surface and a sphere of radius r, is 

2 *^V^-^( P o+ p i£ + ... + P n ? i + ..)df'dt*'. 
Integrating by parts we have 

Therefore the (n + l)t'h term of V is 

- 2 WV»(^+(«-i)^+(«-i)«tt+{»-i)^+i)^+..>^. 



Now R = | ^ + (ij - i-) tfi j ' = a ( 1 + e* p}*) * = a ( 1 - £ /**), rejecting e* 



# 6 2 /*/ 3 



and higher powers of e ; and <p {m) R = <p (m) a — fy^ a . — ^~ to the same degree of 
accuracy. The last member of the expression for V need only be calculated when 
n = ; for all the rest of the terms (^involving / _ l V n dp I where n > o) vanish. 

The first member need only be calculated when n = 0, and when n = 2 ; for when 
w is odd, it vanishes as before ; and also when n is even and greater than 2 : for the 
functions of R involve no higher powers of y) than the square ; and consequently 
they vanish, when multiplied by P 4 , P 6 , . . &c, and integrated with respect to p 1 , from 
— 1 to + 1*. 

When n = 0, the term is 

2 * %f~ i ( R ?/ R - <Pu R ) d V ! ~ 2 ^?/i (r<p,r- <p„r) dp', 
= 4 st g> (a <p, a — <p /y a — -jr a 2 <p a — (r <£, r — (j) n r)L 

* See Art. 9. 
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When n = 2, the term is 

, f ^-j.i/! 1 (^-H(¥+*i?+»V+«-»*^+-)^ 



9> 

2 



9^ 

2 




= 4*^-y)« 8 ^[m-l](^.»-% i f) = -|.» f r«(^-i)««^«. 

To calculate V for the sphere whose radius is r. 
Let the sphere be referred to polar coordinates, the centre 
being the pole. 

LetOC = r,OP=:r 1 ,andPOC=0; thenPC = s/( r * + r* 

Mass of the element at P = g r x 2 <p r Y d r Y sin 9 d 6 d a, and 
Pr p* p2ir prf cj> r x rf^ sin fl dQ d a) __ /^ v^TT ^ ft ri sin 9 dr x d& 

= 2v g J o -^{(r + rj- (r-r l ))dr l = - 7 £ J o r x *pr x dr x 

-ilP(r 2 <p l r-2r<p ll r + 2<p lll r-K), 

K being the value of r x 2 ft r x — 2 rj ft, r x ■+ 2 ft,, r x , when r x = 0. 
The whole value of V then is 

4 r g fafra- </>„«- -i-e 2 a 2 </>« - (<£„r - -f-<fr„r + — ) - ^ r 2 (p* - 3-) e 2 </> a j . 
And the attraction toward the centre 

= -^= 4 ^{^+Mf* 2 -3) re2< H}' 

where 

pr 

$ r = 2 <£„, r — 2 r ^) y/ r + r 2 ^ r — K ==y o r 2 (j>rdr. 

13. To find the attraction of the same spheroid on a particle without it. 
The series (Art. 3.) is 

Now 

A) r > # r /»+ 1 m rf r t _ r /« + l ^ r > _ ( w _j_ ]) r ' w <^ r ' + ^ ( w + 1) fJn - 1 ^ 7 J _ & c ^) ; 

and the general (rath) term of V is 



2 
r 



J^/LV.-i ( r '" +1 ft ^ - (» + r " *« ^ + " (" + ') """ ~ ' +«/ - *"•) rf c' 
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where C is the value of (A) when r 7 = 0. As before, the last term need be calculated 
only when n = 1 ; and the first when n = 1, and n = 3. 
When n = 1, it is 



^/\( R2 ^ R - 2R ^ 



9^ 

2 



When ^ = 3, it is 



+ 24 ^ R) = - y (^ - I) -^ 41a. A 



The whole value of V is 

Ait p 



And 



dV 4* P C f e 2 3 3 / 1W. 1 



14. Instead of the eccentricity e = y/ ( 1 — ^ j it will be more convenient to em- 

ploy the ellipticity e = 1 # . These give s = y And the values of V become for 

an internal point, 

for an external point, 

15. To find the attraction on the supposition that the body is composed of sphe- 
roidal layers, homogeneous in themselves, but differing from one another in density 
and ellipticity. 

First, on an internal point. 

Let /, as before, be the radius vector of any layer ; a 1 its equatorial radius ; g <j> a 1 
its density, and i its ellipticity, being some function of d as % a'. Then 

a! — r> (i + g ' ^'2) an d <£ «' = ^ + r' $' r 1 % r>' 2 = <j> r' + F r> . ^' 2 , suppose. 

Consequently to the term before produced in (12.) by <p r* we must add a term 
similarly produced by F r' . jt*/ 2 . Also, instead of taking, in the first instance, the por- 
tion comprised between the surface and a sphere of radius r, we must take the por- 

* See Puissant, vol. i. p. 259 ? where the word ellipticity is used in this sense. 
MDCCCXLI. N 
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tion between the surface and that spheroid on which the point lies, whose ellipticity is s x . 
The (n + l)th term of V now becomes 

- (a <p i a - <p u a ~~ ^ a 2 <p a - ~|(^ 2 — j) r 2 <p a) 



The first part of this gives 4wg 



a being semiaxis major of the stratum on which r lies. To determine the other part 5 
it is necessary to compute it when n = and n = 2, which gives 

To the sum of these we must add the value of V for the inner spheroid ; and for this 
purpose we have to obtain V for an external point. 
To the expression in (13.) we must add 

2 

r 



2 -^f\dp< (^ P„_ J^ FrV» + * rfr»), 



to be calculated when n = 1 and n = 3. This is 



4t? 



{rO" F »' . 



a /2 ^ f + 



F 3 - 



5r 8 



fFaf.tfrdafl. 



The whole value of V is 



4f 



^4^a4^i)^ a +?r Fau ^ a ' 



^ 



2 _ 



>. 



+ 



f a ¥a! .a^da! 



.J 



After writing a for a, and s x for s, add this to the other value of V 5 and apply the equa- 
tions 

7-jC" F d • a ' 2 <* a ' = t£ a ' 3 id $ d = a2 h P a - ±f o * pa' d (tfi *'), 
and similar equations for the other integrals ; and we shall obtain 



V = 4Tf 



a^a-^a-a^a+p^a+^-^-^^tpa'dia'h') +^f \a'd{a!h')\ 
~ T (f 2 " t) {£'<"*<" + ?/***<* d(*i) } 
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16. To find the equation of equilibrium of a heterogeneous spheroidal mass of fluid, 
revolving about its axis, with an angular velocity *. 

By the principles of hydrostatics the general equation is J ~t =J (X dx + Y dy 

+ Z dz), § r being the density at the point (x,jj 9 z), p the pressure, and X, Y, Z the 

dY dY dY 

sums of the resolved parts of the forces ; which are — -^, — ^--, — -^ and the 

centrifugal force. Let the axis of z be that of rotation ; then the centrifugal force is 
m % x along x, and ® 2 y along y. Let us express ® in terms of the ratio of the centri- 
fugal force at the equator to the equatorial gravity. Call this ratio m, which is small 
in the case of the earth, being of the same order as g. Then 



Therefore 



o? a 3 w a s 4: worn fy a 

mass 4 w g ty a p a 8 



Y dY 4vgmtya.x v dY 4w§ m^a.y ^ • 



dV 



rf^ T a 3 * dy ~ a 3 * dz> 

and 

Now^/ ~~f is a constant for a level surface. Hence for any stratum we have 

c = V + £Ef£±f (1 _ ^) ,a. 
At the surface this is 

C = *? - T^O " (* - f 0* - T>?/' + * * ^ + 1 ^? (l - ^V, 

where 

/* a . i ^ (a' 8 (1 - * f )) ir i xt r a ± t d ( a *l) i i 

Ma=zJ o tfx*'- 1 -^? ^da f ,andNa = y o $ a 1 -~^r~ f d a 1 , 

For r write a (1 — g f£ 2 ) 3 then 

^ Ma ,_ . 9X 1 / 9 l\Na 1 , 9N Ma 

Equate the coefficients of ft 2 , then 

Ma/ m\ 3Na /TIN 

~T V 1 - TV = -53r----( B )- ' 

17, By differentiating and changing the sign of j-jr, we obtain the amount of 

gravity which acts towards the centre ; which, to the order we are now considering, 
is the same as the whole force of gravity ; since the cosine of the angle of the verti- 
cal differs from unity only by terms of a higher order. 

n2 
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Consequently 

("Ma %( 9 l\Na m , .rMal 

CM a „ . ■ ox 3 / „ l\Na m ,_ > Mai 

= 4» f | 5p -(l + 2i^- T (^- T j- 3 r- T (l-^)- 5 rj, 

= 4* f ^{l + 2.^ + (l-3^)(«-*)-«(l-^}(byB.), 

= 4^^{l+ g - 3 -f+^(|m- s )} = G{l +S in 2 /(4 W - S )}, 

where G is the equatorial gravity and / the latitude. 

18. From this it appears that up to terms of the 1st order, R = a (1 — i sin 2 /) is the 
equation of the curve which generates the surface of equilibrium, where the value of 
s depends on m, or on the velocity of rotation : but as the coefficients of higher 
powers of sin I may be considerable, it will be useful to find the surface of equilibrium 
to a greater degree of exactness. For this purpose we must introduce the fourth 
power of sin l^ whose coefficient will be of the second order. Let the equation of the 
strata be d = r' (1 + s> 2 + A 1 ^), or r f = d (l - s> 2 + (s f2 — A') ^ 4 ), where s' and A' 
are functions of d as % d and d. Then 

^ a f — <|> r » + p r ^2 + r 2^Clli£r) £ + r ^ r 0r^ 4 = ^r+Fr.^ 2 + nr.^ 4 ? suppose. 
The (V* + l)th term of V for an internal point now becomes 

By writing for r and R their values in terms of a and a, it is easily found that 
X* J& d r- equals ^ « + (-£,)(- . ,. + ^. ( , _ A , _ p , !=• ,,) 

+ ( n -3") g 2 — the same functions of a and s x .- 

And that similar equations are true for the two remaining integrals. 
1st. Let n = 3 and we have 



*gj-i d t*'' 



f"<£ 4> a' d a! + a 2 </> a ( - s p' 2 + ^ (-| s^- A)) + a 3 </>' a-%- 

+ p' 2 ^" a' . F a' rfa' + p* a 2 F a (- c /»'») + i *' 4 % /*V n a' da' 

. — same function of a and s l3 
which gives 

4 *f $J^et<pddei + a 2 pa (- y + y(-§-s 2 - A)) + a 3 p'a^ + \ J^Fa! .a! .dcf 



— « 8 F a y 4- yy a' . II a' rf a' > — 4 ff f (same function of a x gj and Aj). 



AND THE FIGURE OF THE EARTH. 



93 



6 

~2 



n 



2nd. Let n — 2, and we get 

-f {a ?' a - z<p' z)*^ +fj^ d a! + {F a - F &) {- J p'*) +J^Tld .a'.d 



A')) 



a 



which gives 

4ir e (^-y)^a(-y+^(4-A)) + fl^a^s2+^ ^fdd 



«Fa' 



6 



35 



;Faj 



6 /•* n d I / 1 \ 

4- or / —j-dd \ — 4 7T g> I ^ 2 —■3"/ r2 ( same function of a, s x and A x ). 



3rd. Let 72 = 3, and we have 
1225 



64 



(/ 



6 

7 



p s +IK^/-\^->h4)U"^«'+(¥ 



$ a 
a 



) 



(-^»- f « ( ? + .)) + (^4)4 + /'5« + ( 



Fa Fa 



a 



2 



a 



2 



)( 



g' jh/ 2 ) 



+y a 



«IIa' 



a' 



/3 



^a' 



} 



or 



( 



6 



f* 



4_-^,,2_|_^ 4 ^ ?r 4 



7 



f*' 



<f> a 



a 



2 



t(t + a) + *i 



Fa 



8 9 a 



2 



+ 9 Jo a® d 



a 



-( 



pfi 



1 

7 



3 \ 
[jb 2 + gg ) 4 w § r 4 (same function of a, s } and A a ). 



19. To the sum of these must be added V for the inner spheroid^ for which we shall 
have to find the value of V generally, when the attracted point is without the body. 
The general term of the series for V is 

?±Jl x dfb' P n _ A J o r" 8 .<pr' .dr' + (i!*J o r ,n Fr< . dr 1 + ^J q ^ .Ur' .dr' i , 



r 

and 



J %R r ,n <pr ! dr ! = fj d n .<pd .dd + a*+ l <pa (- iy? + ^ (s 2 - A) + ^4 -J g«) 



and similar equations hold for the other functions. 
1st. Let rc = 2, then we have 

2 



~f-\ dp 1 \ Jl* d 2 <pd dd + a 3 <p a (- g ^ 2 + ^ (2 s 2 - A)) + a*^ 



.2 



a 



10 



which is equal to 

^5 f Q \ ,2 <pddd + a*pa(---l~ + y(2g2-A))+a*?/fl^+ ^-J^'aVFa'.d 

— y a 3 F a . s + -jr J d 2 Tld dd > . 



a' 
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2nd. Let n = 4 3 and the term is 



!9 

4 



+ cfi a^f +fj ci* ¥c£dct-a*Ya* ^ + fj d* n a 1 da! j , 



or 



4ng 



~§~ 5 aSFa.s+fJ o a amaida'}« 
3rd. Let n = 6, and we get 

i 5 o* - 1 ^ + 1) 4 -?/> "'("'*- f ^ + a) { r ^ ^ " 

+ d?<pa(~i^ + p'4 (4 a 2 - A)) + a* <p' a ^|- +J^'d 6 Fddd -a?Fas pF 

+f o a a'nia<da!'j, 
which is 

20. In these three terras write a for a, i x for e, and Aj for A, and add them to the 
other value of V, and apply the equations 

fjd* n d dd = yX" a ' 4 d (X «' 2 ?' a ') +X" a ' 3 ^«' ^ (^ «') = y « 4 j2 ?>' " + a 3 A ^ a 
— 2f Q a d s %di.d<pd —f" <pdd (d 3 A') = y a* s 2 p' « + a 3 A p a — 2 a 3 s 2 © a 

+ 2^" p a' d (a* s' 2 ) -fjtpdd (a' 3 A'), 

and similar equations for the other integrals, and the value of V for an internal point 
becomes 

y a ^ — L dd — ± - da '+s-Jo * a — 1 —m J ~ da 

-5^ 2 -F) V r Va *<* dd da +^Jo *<* - L " rf3 ~ ^7 



h 



+ 9V' 4 ~ J^+ 35) V 4 Ja *<*— -33 M+7J ^'- da' 



d 




- 4w S $M'a- M'a + ^Ma + y (f* 2 - T ) (r 2 N'« - r 2 N'a + 75 Na) 

+ i(f^ 4 - y ^ 2 + i) (r*Pa - r*Pa+^Pa) }, 
where the functions M, N, P> &c. are substituted for the corresponding integrals. 
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21 • The equation of equilibrium will be formed exactly as In (16.), except that 
the expression for the velocity must be found more accurately. By (17*)> the force 

of gravity at the equator is 4 te § ^ ( 1 + s 5- ) ; and the centrifugal force is a> 2 a ; 

therefore m = — - WT ~, ^-, and «« = 4 * g ^ (w + m * - -g-J. 

The equation of equilibrium becomes 

/£ = c-v + », f ,.(i _„£„(,+,_»•), 

which at the surface becomes 

_,, Ma 1 / „ 1\ N« 1 / 6 ,. 8\P«. 1., 9N ,M«/ , 3w 2 \ 

+ ^_„ ( ! _ 2 g ^2) ^ 3 ^ m + m § _ _j . 

Equate the coefficients of ^ 2 and those of ^ 4 5 and we have two equations for deter- 
mining s and A, thus showing* that equilibrium is possible. These are. 

Ma/ m 3 3 A Nfl/1 s\ 2 Pa 

_ ^ ._ _ _ _ OTi + _^2j _ _. (_ _ _) _ __ = 05 

, Ma, A Na 3s , 1 Pa 

and _(A + i W i)- ? . T + - 7 .=0; 

Mb/ m , 6A , 3 , ms 11 „\ Na 

whence _^_ T+ _ + _ m 2 + _. yi 2j == _ ?f 

Ma / 5ros\ 1 Pa 

22, The resultant attraction in the direction of r, is obtained by differentiating 
J—r, as found in (21.), and changing its sign. This produces 

A fMa 3/ l\Na.5/ . 6 9| 3\Pa . „. / f 3rf\Mal 

since the terms arising from differentiating the expression for V in (20.) with respect 
to a 5 cancel each other. 

For r write its value, and arrange the result according to powers of \h, and this 
formula becomes 

Ma/. . 3 m 2 

3 



- 2 ( 1 — m + m i jj- +f 2 (2s + ffl-f2«»s a~) +j!4 4 (2A + j 2 - me) ) 

N« / 1 2/ 3 4 x .12\ . Pfl/1 2 10 . * 5 \ 

~ IF I" 5 +f* (t - TO +^ 4 sT 8 ) + -^ V21 ~ P* 2 2l + f* 4 9J 

' ' Maf 3m 27 , 9 2 , 3 

+ ^ (|. m _ g + 7_ 2 w g _ ^ m2 _ ^ g2 + ^ A ) + ^4 ( 4 g2 _ 3 A _ ^ m g ) I ^ 



96 MR. C. J. HARGREAVE ON THE CALCULATION OF ATTRACTIONS, 

by applying the two equations deduced in last article. The total attraction normal 
to the surface, or the force of gravity, is found from this by dividing by the cosine of 
the angle of the vertical, or by 1 — 2 s 2 {(jfi — p 1 ) ; whence 



M a 



1 + g g jgm e + ^ m 2 — j s 2 + j A 



Let X be the sine of the real latitude ; then as (i 2 = X 2 — 4 e X 2 — X 4 , we get 



g = 4v S 



Ma 



3d 



2 



1. 3m ^ . 9 2 2 2 . 3 A 

l_f_ g 2 j^me + ^m 2 — -yg 2 + yA 



+ X 2 (4m- s + ym g -^m 2 +ye 2 + yA)-X4(2 S 2 +3A-4m g )]. 
Let G represent the equatorial gravity ; then 
£ = G{l + X 2 (f m- s-y me + y S 2 + yA)~ X*(2e 2 + 3A-f- ms)j 

= GJl + sin 2 z(|m-s- Y^mi + y g 2 - y a) + sin 2 Zcos 2 1(2 s 2 + 3 A- |to«) j, 

which is an extension of Clairaut's theorem. 

23. In this process A indicates the amount of deviation of the required surface 
from the surface represented by r = a (1 — s ^ 2 ). If the equation had been taken 

a = r < 1 + e p 2 — (^ ^ e 2 + BJ f f^ 4 — ^ 2 j -f e 2 f£ 4 > , B would have been the index of 

deviation from an elliptic spheroid. 

To apply Clairaut's theorem to this surface, we have 

3 e 2 

g = e + g- e 2 + B, and A = — ^ — B ; 

whence 

g = G jl + sin 2 / (■§- m -e + jB-^me) - sin 2 /cos 2 / (| me - ^ + 3 a) j, 
which is the same expression as that obtained by Mr. Airy in the Philosophical Trans- 
actions, 1826 5 except that instead of e or — — , the symbol is used to represent — ^- . 

24. The circumstance of the terms arising from the differentiation of V with respect 
to a vanishing, affords an easy method of extending Clairaut's theorem indefi- 
nitely, without calculating the value of V. 

It may be shown independently, that these terms cancel each other in all cases. 
The (w+ l)th term of V for the portion including the point is 



The corresponding term for the other portion is 

2 
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and it is evident, by inspection of these functions, that that portion of the sum of their 
partial differential coefficients, which arises from differentiating with respect to r 
under the signs of integration, is equal to zero. This being the case, it is not neces- 
sary to know the forms of the functions M a 9 N a, P a, &c, nor their numerical coeffi- 
cients; but only the functions of /a, by which they are respectively multiplied; and 
these are Laplace's coefficients. 

Thus the equation of equilibrium at the surface would be 

Ma 1 / 9 1\N« , 1 / . 6 , 3\Pa / , 15 . , 5 _ 5 \Qa 

l ii 

+ 2~ (1 — ^ 2 ) r 2 m G, where — = — (1 + s ^ 2 + A ^ + D ^ 6 ), suppose. 

Expand r, recollecting that N a, P a, and Q a are of the 1st, 2nd, and 3rd orders 
respectively, and we have three equations to determine s, A and D. 

By differentiating C with respect to r, and eliminating N a, P a, and Q a by these 
three equations, we have the resolved force in r, which divided by the cosine of the 
angle of the vertical gives g exactly as in (23.). 

It is evident that this may be carried on indefinitely ; and to any order, without 
finding g for the next lower order. 
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